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BASIC  ANALYTICAL  MODEL  FOR 
ENVIRONMENTAL  IMPACT  ASSESSMENT 
OF  SURFACE  WATER  RESOURCES- 
DOSAG  USER  MANUAL 


1 INTRODUCTION 


Background 

The  Environmental  Impact  Computer  System 
(Lies),'*  developed  by  the  U.S.  Army  Construction 
Engineering  Research  Laboratory  (CERL)  provides 
systematic  procedures  and  a technical  data  base  to  help 
Army  planners  and  decision-makers  perform  environ- 
mental analysis  and  prepare  and  review  environmental 
impact  assessments  and  statements  (EIAs/ElSs).  EICS 
provides  qualitative  information  but  does  not  provide 
the  tools  needed  to  quantify  environmental  impacts. 
The  Council  on  Environmental  Quality  requires  quanti- 
fication whenever  possible  for  ElAs/EISs. 

The  computerized  model  discussed  in  this  report 
was  developed  to  help  fulfill  the  need  for  quantitative 
water  quality  analysis  during  environmental  impact 
analysis.  It  will  be  used  primarily  by  Army  planners 
and  decision-makers  associated  with  Facility  Engineer- 
ing to  obtain  information  for  environmental  analysis 
preparatory  to  writing  an  EIA/EIS. 

Objective 

The  objective  of  this  study  was  to  develop  and  pro- 
vide user  instructions  for  a computerized  basic  water 
quality  model  which  could  be  used  to  quantify  impacts 
on  the  dissolved  oxygen  content  of  a stream  resulting 
from  Army  activities.  This  model  can  be  used  for  a 
single-reach,**  single-point-source  pollutant  inflow 
with  steady-flow  conditions. 

' L.  V.  Urban,  et  aL,  Computer- Aided  Environmental  Impact 
Analysis  for  Construction  Activities-  User  Manual,  Technical 
Report  K-50/ADA008988  (U.S.  Army  Construction  Engineer- 
ing Research  Laboratory,  March  1975). 

•EICS  helps  Army  planners  efficiently  identify  primary  and 
secondary  impacts  of  proposed  Army  projects  or  activities  and 
suggests  ways  to  mitigate  these  impacts.  EICS  is  part  of  the 
Environmental  Technical  Information  System,  which  includes 
two  other  operating  systems:  the  Economic  Impact  Forecast 
System  and  the  Computer-Aided  Environmental  Legislative 
Data  System. 

••Section  of  a stream  unaffected  by  water  withdrawals 
(water  treatment  plant),  discharges  (waste  treatment  plant 
effluents,  tributaries),  or  other  physical  modifications  like 
dams,  weirs,  locks,  etc. 


Approach 

Three  analysis  techniques  are  combined  into  a single 
model  capable  of  quantifying  impacts  relating  to  DO, 
and  related  parameters  of  BOD  and  water  temperature. 
Fourier  series  analysis  is  used  with  historical  water 
temperature  data  to  develop  an  annual  water  tempera- 
ture profile.  Five-day  BOD  data  are  used  to  calculate 
the  ultimate  BOD  and  deoxygenation  rate  coefficients 
for  both  carbonaceous  and  nitrogenous  waste  discharges 
for  receiving  water.  The  effects  of  reaeration  and  photo- 
synthesis are  analyzed  in  the  computations. 

Model  Requirements 

Several  factors  were  considered  when  selecting  tech- 
niques for  predicting  the  effect  of  waste  effluents  on 
the  dissolved  oxygen  (DO)  content  of  streams.  The 
modeling  of  DO  in  streams  is  a well-developed  science, 
and  several  models  can  be  used  for  this  purpose. 
DOSAG-1  and  QUAL 1 are  examples  of  well-documented 
models  which  have  been  proven  througli  widespread 
application;^  however,  these  models  are  more  complex 
than  appears  to  be  necessary.  Two  or  three  weeks  of 
data  acquisition  and  preparation  are  often  needed  to 
use  these  types  of  models.  Clearly,  Army  installation 
personnel  would  have  difficulty  conducting  such 
analysis  on  a routine  basis. 

A variety  of  situations  can  result  in  stream  BOD 
loading.  The  most  familiar  is  point  source  effluent 
from  wastewater  treatment  plants.  Area  source  pollu- 
tion has  received  great  interest  recently;  however, 
this  type  of  analysis  is  also  very  complex.  For  impact 
analysis,  techniques  which  consider  the  primary  DO 
sources  (reaeration  and  photosynthesis)  and  the 
primary  DO  sinks  (carbonaceous  and  nitrogenous 
biochemical  oxygen  demand  [BOD])  are  adequate. 
The  usual  situations  for  which  these  techniques  will 
be  used  and  the  capabilities  of  the  typical  user  sug- 
gest that  a single-reach,  steady-flow,  point  source 
model  is  appropriate.f  The  techniques  provide  rapid 
analysis  for  users  having  little  expertise  in  water- 
quality  modeling  and  require  a minimum  of  input  data. 
The  most  common  use  of  these  techniques  is  expected 
to  be  prediction  of  impacts  from  various  levels  of 
wastewater  treatment  plant  operating  efficiency.  The 
worst-case  situation  is  expected  to  occur  during  low- 
flow,  high-temperature  stream  conditions. 

^Simulation  of  h'ater  Quality  in  Streams  and  Canals.  Pro- 
gram  Documentation  and  User  Manual  (Texas  Water  Develop- 
ment Board.  September  1970). 

tVerification  of  this  approach  is  under  way.  Techniques  to 
expand  the  model  to  include  additional  tributaries  and  effluent 
sources  are  being  investigated. 
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Mode  of  Technology  Transfer 

User-related  infurniation  from  this  report  will  be 
issued  as  an  Engineering  Technical  Letter.  The  com- 
puter program  will  be  documented  in  another  report, 
and  the  source  deck  will  be  available  from  National 
Technical  Information  Service.  The  computer  program 
will  remain  operative  on  the  time-sharing  computer 
system  (Naval  Shipyard  Research  and  Development 
Center)  which  would  allow  access  by  field  users. 

After  the  decision  is  made  concerning  where  the 
Environmental  Technical  Information  System  (ETIS) 
will  be  field-implemented,  the  model  will  be  integrated 
and  appropriately  referenced  in  the  ETIS. 


2 ANALYSIS  FUNCTIONS 

DOSAG,  the  computer  program  for  this  model,  is 
a collection  of  subroutines  wltich  perform  basic  water 
quality  analysis  functions,  and  is  used  to  quantify 
impacts  on  BOD  and  DO.  DOSAG  has  been  developed 
to  evaluate: 

1 . Annual  variations  in  a water  body  temperature 
through  Fourier  Series  analysis,  using  historical  data 

2.  The  ultimate  BODs  and  deoxygenation  rate  co- 
efficients for  both  carbonaceous  and  nitrogenous 
waste  discharges 

3.  The  impact  of  carbonaceous  and  nitrogenous 
waste  discharges  on  the  DO  of  a receiving  water  body. 
The  program  has  several  options  for  reaeration  or  photo- 
synthesis. 

Temperature 

Temperature  is  probably  one  of  the  most  significant 
and  widely  measured  parameters  among  the  water 
quality  characteristics  of  surface  waters.  Most  of  the 
physical,  chemical,  and  biological  properties  of  water 
are  functions  of  temperature.  As  temperature  increases, 
the  saturation  concentration  of  DO  decreases  and  the 
atmospheric  reaeration  rate  increases.  An  increase  in 
water  temperature  has  both  beneficial  and  deleterious 
effects  on  the  environment.  Studies  have  indicated  that 
for  each  IO°F  increase  in  temperature,*  the  cost  of 
chemicals  to  coagulate  and  treat  water  effectively  can 


*0.  E.  Arnold,  ‘Thermal  Pollution  of  Surface  Supplies,” 
Journal  American  Water  Works  Association.  .S4 : 1 332  (Novem- 
ber 19621. 


be  reduced  by  $0.30  to  $0.50  per  million  gallons.  The 
effectiveness  of  chlorination  and  the  efficiency  of 
sedimentation  tanks  and  filtration  units  are  enhanced 
by  higher  temperatures.  On  the  other  hand,  high  tem- 
peratures increase  metabolic  rates  of  organisms,  which 
increases  the  demand  for  DO.  Higher  temperatures  may 
also  promote  algae  growth  if  other  conditions  are 
favorable. 

Traditionally,  analytical  models  for  water  tempera- 
tures have  been  based  on  two  different  precepts.  The 
first  is  the  heat  budget  of  a body  of  water,  which  con- 
siders factors  like  incident  solar  radiation,  atmospheric 
radiation,  intensity  of  reflected  solar  radiation,  back 
radiation  or  radiation  emitted  from  the  surface  of  the 
water  body,  conducted  heat,  loss  of  heat  due  to  evapo- 
ration, and  convection  of  heat  with  evaporated  water, 
etc.  This  method  requires  extensive  instrumentation  to 
evaluate  various  factors  affecting  the  heat  budget  and 
is  prohibitively  expensive  and  time-consuming. 

The  second  approach  is  to  develop  an  analytical 
model  for  water  temperature  based  on  historical  water 
temperature  data  for  any  stream  or  lake.  The  existence 
of  cyclic  variations  in  mean  daily,  weekly,  or  monthly 
water  temperatures  was  reported  by  Kothandaraman 
and  Evans,^  Thomann,*  and  Ward.®  These  invesxigators 
indicated  that  the  harmonic  function  with  a periodicity 
of  1 year  fits  water  temperature  data  in  lakes,  estuaries, 
and  rivers,  and  accounts  for  approximately  95  percent 
of  the  total  variance. 

The  annual  temperature  variations  have  been  found 
to  fit  an  empirical  equation  of  the  form: 

T = p sinbx  -I-  q cosbx  + T [Eq  1 1 

where  T is  the  calculated  water  temperature  (°C),  x is 
the  number  of  days  since  January  l(x=  Ifor  January  1), 
p and  q are  the  components  of  the  harmonic,  b is  a 
constant  given  as 


^V.  Kothandaraman  and  R.  L.  Evans, ‘‘Annual  Temperature 
Variations  in  an  Impoundment  in  Central  Illinois,”  Journal 
American  Water  Works  Association,  62:639  (October  1970). 

*R.  V.  Thomann,  “Time  Series  Analysis  of  Water  Quality 
Data,”  Journal  of  the  Sanitary  Engineering  Division.  93:1 
(American  Society  of  Civil  Engineers  [ASCE],  February 
1967). 

®J.  C.  Ward,  “Annual  Variation  of  Stream  Water  Tem- 
perature,” Journal  of  the  Sanitary  Engineering  Division. 
89:1  (ASCE,  December  1963). 
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b = days  = 0.987  degrees  per  day  or  [Eq  2) 
0.0 1 72  radians  per  day 

and  T is  the  mean  annual  temperature  of  the  observed 
record. 

The  standard  deviation,  a,  for  a set  of  observations 
is  given  by  the  following  equation; 


The  standard  error  of  estimate,  S , is  given  by 


/S(Ti-T)^ 


and  the  correlation  coefficient,  r,  is  given  by 


where  Tj  is  the  observed  temperature  and  N is  the  total 
number  of  observations.  N pairs  of  evenly  spaced  val- 
ues of  X and  Tj  covering  the  entire  period  of  1 year  of 
multiples  thereof  are  required  to  use  these  algorithms. 
Figure  1 shows  various  entities  of  this  temperature 
model. 


The  parameters  p and  q of  Eq  1 are  determined  so 
that  the  sum  of  the  squares  of  the  residual  errors,  , 
i.e.,  the  sum  of  the  squares  of  the  differences  between 
the  predicted  and  observed  temperatures,  is  a minimum. 
This  is  expressed  mathematically  as 

= 2 (p  sinbx  + q cosbx  -f  T — Tjl  ^ [Eq  6] 

This  condition  is  satisfied  if  the  first  partial  derivatives 
of  R with  respect  to  p and  q are  zero. 

^/yD2\  _ 

— — = 22(p  sinbx  + q cosbx  + T — Tj)sinbx  = 0 

[Eq7] 

— — = 22(p  sinbx  + q cosbx  -I-  T — Tj)  cosbx  = 0 

(Eq  8] 

Eqs  7 and  8 can  be  rearranged  to  give  the  two  normal 
equations  for  evaluating  p and  q as  follows; 

pSsin^bx  + qSsinbx  cosbx 

+ T Isinbx  — I TjSinbx  = 0 [Eq  9| 


pSsinbx  cosbx  4-  qScos^bx 

+ T Scosbx  —2  TjCOsbx  = 0 [Eq  1 0] 


t 

K. 

11“ 

Ji!  .-L  ^ 


^ r+Sy 


Ji4N  F£B  MAR  APR  MACf  JUN  JUL  AUG  SEP  OCT  NOV  DEC  JAN 


Figure  I.  Annual  variations  in  water  temperatures. 
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Solving  these  two  equations  for  p and  q provides 

^ Idunbx  iosbiiHi;T|SinbJi  Tlsinbx)  isin'bxCS;  l.cosbx  Ticosbx)! 
^ (isinbx  cmbx)^  — (Isin^bxMXcos^'bx) 

[Eq  1 1 1 
and 

^ (2T,sjnbx— Tiisinbx)  Icos^bx  (ITiCosbx-f  icosbx)isinbx  cosbx 
^ i^sin^bx  icos^bx--(lsmbx  cosbx)^ 

(Eq  12] 

Biochemical  Oxygen  Demand 

The  wastes  discliarged  into  a stream  from  municipal 
and  industrial  plants  contain  a variety  of  chemical  com- 
pounds. Of  primary  interest  is  that  portion  of  the  waste 
which  is  biodegradable,  and  hence  oxygen-consuming. 
When  this  material  is  introduced  into  a water  course,  it 
undergoes  biochemical  oxidation  when  micro-organisms 
use  the  organic  matter  for  energy  and  growth. 

If  enough  DO  is  present  in  the  water,  the  microflora 
will  be  primarily  aerobic,  using  the  DO  to  carry  out 
oxidation  reactions,  and  producing  water  and  carbon 
dioxide  as  end  products.  However,  if  enough  oxygen  is 
not  present,  anaerobic  organisms  predominate,  provid- 


ing undesirable  end  products.  In  addition,  the  concen- 
tration levels  of  oxygen  have  profound  effects  on  the 
physiology  and  type  of  fish  found  in  the  water  bodies. 
The  lethal  effect  of  low  concentrations  of  DO  appears 
to  be  increased  by  the  presence  of  toxic  substances, 
such  as  excessive  dissolved  carbon  dioxide,  ammonia, 
lead,  copper,  etc. 

Because  of  the  variety  of  oxygen-demanding  organics 
in  the  waste,  it  is  common  to  measure  the  strength  of 
wastes  in  terms  of  BOD  rather  than  to  analyze  the 
chemical  constituents  of  the  wastes.  Since  the  BOD 
reaction  is  closely  related  to  the  unimolecuiar  type  of 
reaction,  a plot  of  the  amount  of  the  organic  matter 
remaining  versus  time  yields  a parabolic  curve  similar 
to  the  decay  curve  for  a radioactive  element  (see  Figure 
2).  Likewise,  a plot  showing  the  amount  of  organic 
matter  oxidized  versus  time  produces  a parabolic  curve 
that  is  the  reciprocal  of  the  first. 

Since  the  amount  of  oxygen  used  is  directly  pro- 
portional to  the  amount  of  organic  matter  oxidized  in 
biochemical  oxidations,  a plot  of  oxygen  used  versus 
time  would  produce  a parabolic-type  curve  like  the 
one  produced  for  organic  matter  oxidized  (Figure  2). 


TtME,  DAYS 


Figure  2.  Plots  of  oxygen  used  and  organic  matter  oxidized. 


10 


Figure  3.  Typical  BOD  curve. 


Figure  3 is  a typical  BOD  curve  fora  municipal  effluent. 
For  the  first  8 to  10  days,  its  characteristics  are  similar 
to  those  of  the  :jrve  for  organic  matter  oxidized 
(Figure  2).  After  that  time,  the  BOD  curve  digresses 
radically  from  the  course  it  would  be  expected  to  follow 
as  a monomolecular  or  first-order  reaction.  This  is 
because  the  nitrifying  bacterial  population  has  become 
large  enougli  to  significantly  increase  the  oxygen  de- 
mand exerted  by  those  bacteria.  These  organisms 
oxidize  nitrogen  (in  the  form  of  ammonia)  to  nitrite 
and  nitrate. 

On  the  basis  of  first-order  kinetics,  the  rate  of  BOD 
removal  in  a water  body  is  given  by 

^=-K,L,-K„L„  [Eqi31 

where 

L = total  BOD  remaining  at  any  time,  t,  mg/? 

Lj,  = carbonaceous  BOD  remaining  at  any  time, 
t,  mg/8 


L„  = nitrogenous  BOD  remaining  at  any  time,  t, 
mg/8 

Kg  = carbonaceous  deoxygenation  coefficient  to 
the  base  e,  per  day 

K„  = nitrogenous  deoxygenation  coefficient  to  the 
base  e,  per  day 

Integration  of  the  carbonaceous  term  of  this  equation 
yields; 

Lc  = L„,e-'^'  [Eql4] 

where 

L„g  is  the  ultimate  carbonaceous  demand.  Integration 
of  the  nitrogenous  term  of  Eq  13  yields: 

L„  = Lure”""'  [EqlS] 

where  is  the  ultimate  nitrogenous  demand. 

If  L=  Eg  -F  L„,  then 

1 = L e“'^c'  + I 

^ i-uc  •'un 


n 


lEq  16) 


If  there  is  a time  lag,  tg,  for  the  onset  of  nitrogenous 
demand  from  the  carbonaceous  demand,  which  invari- 
ably is  the  case,  then  Eq  16  takes  the  form 

L=L„,e-''c«  + (Eq  17| 

In  the  above  equation,  the  factors  Lyp,  K^,  and  Kp 
must  be  determined  from  the  BOD  time  progression 
curve  obtained  under  standardized  laboratory  condi- 
tions. To  differentiate  between  the  carbonaceous 
oxygen  demand  and  the  nitrogenous  oxygen  demand, 
the  laboratory  conducts  tests  on  waste  stream  samples 
both  with  and  without  adding  inhibitory  agents  to 
suppress  the  activity  of  nitrifying  bacteria.  BOD  pro- 
gression curves  obtained  from  the  samples  without  an 
inhibitor  yield  the  combined  BOD  resulting  from  the 
oxidation  of  carbonaceous  matter  and  ammonia. 
Results  from  samples  containing  an  inhibitor  yield  the 
BOD  values  resulting  from  the  oxidation  of  carbona- 
ceous matter  alone.  The  difference  in  the  values  obtained 
under  these  two  conditions  yields  information  on  the 
BOD  resulting  from  nitrification.  Various  chemicals  have 
been  tried  as  inhibitors;  methylene  blue,  allylthiourea, 
and  acidification  are  among  the  most  widely  used. 

Ultimate  oxygen  demands  and  deoxygenation  co- 
efficients are  estimated  from  the  BOD  progession  data 
by  means  of  analytical  procedures  The  Thomas  slope 
method,  the  method  of  moments,  the  rapid  ratio,  the 
daily  difference  method,  the  Reed-Theriault  method, 
and  the  least  squares  method  of  nonlinear  curve  fitting 
have  all  been  used  to  estimate  ultimate  BOD  and  the 
deoxygenation  coefficients.  Each  method  gives  a slightly 
different  set  of  values  for  these  two  parameters.  The 
Reed-Theriault  method,  which  provides  reliable  infor- 
mation of  ultimate  BOD  and  deoxygenation  rate  con- 
sistent with  the  effort  involved,  is  explained  in  the 
following  paragraphs. 


is  applicable  equally  to  data  for  both  carbonaceous 
demand  and  for  nitrogenous  demand.)  This  procedure 
is  then  used  together  with  the  observed  data  to  deter- 
mine a correction  factor  to  the  assumed  K value.  The 
operation  is  repeated  with  the  corrected  K value  until 
the  correction  factor  is  small,  i.e.,  until  the  assumed 
and  calculated  values  of  K approach  one  another.  The 
details  of  the  Reed-Theriault  algorithm  are  given  in  Eqs 
18  through  26. 

Y = Lp(l-e-'")  (Eql81 

where  Y is  the  BOD  satisfied  in  time,  t,  and  all  the 
other  terms  are  as  defined  earlier.  For  the  purpose  of 
reducing  this  expression  to  a linear  form,  an  approx- 
imate value  of  K has  been  assigned  so  that  the  most 
probable  value  of  K is  given  by  the  relation 

K = K*  -I-  h [Eq  19] 

in  which  K'  is  the  trial  value,  and  h is  a relatively  small 
number,  positive  or  negative,  whose  magnitude  is  to  be 
determined.  Substituting  this  value  into  Eq  18, 

Y=  Lu[l-e-('^''^'’>‘l 
= L„a-e-’''’)+L„hte-'^‘‘  [Eq  20] 

This  is  of  the  form 

Y = af,  -t-bfj  [Eq21| 

in  which  a = L^,  b = Lyh,  fi  = 1— e"*^'*,  and  fj  = 

te"*^  *.  For  Eq  21,  the  difference  between  an  observed 
and  a calculated  quantity  is  Y — afj  — bfj . The  sum,  R, 
of  the  squares  of  the  residual  errors  is 

R=  2(Y-af, -bfj)^  [Eq22] 


The  Reed-Theriault  method"'  estimates  the  ultimate 
BOD  and  deoxygenation  coefficient  from  the  observed 
BOD  progression,  using  the  least  squares  principle  to 
minimize  the  residual  errors  between  the  observed  data 
and  the  fitted  monomolecular  curve.  To  apply  the  pro- 
cedure, a method  of  successive  approximations  involving 
trial  and  error  is  used  in  which  it  is  necessary  to  assume 
a trial  value  of  K.  (Note:  The  subscripts  of  c and  n for 
K used  earlier  have  been  dropped  here.  The  procedure 


The  most  probable  values  of  the  constants  a and  b will 
be  obtained  when  R is  a minimum.  Equations  defining 
this  condition  are: 

|^=2:f,Y-a2:f?-b2f,f2  =0  [Eq  23] 
II  = ZfjY-aSfifi -b2f|  =0  [Eq  24] 


^Lowell  J,  Reed  and  Emery  J.  Theriaull,  “The  Statistical 
Treatment  of  Reaction-Velocity  Data  II,  Least  Squares  Treat- 
ment of  the  Unimolecular  Expression:  Y = L 
Journal  of  Physical  Chemistry,  Vol  35(4)  (1931),  pp  950-971. 


The  solutions  of  Eqs  23  and  24  are: 


2fl  2f,Y-  Sf.fj  SfjY 
Sf?  Zfl  -(2f,f 


|Eq2S] 
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Sf?  SfjY-Sf.fj  2f,Y 
Zf?  2f| -(Sfifi)* 


[Eq  261 


From  the  values  of  a and  b as  obtained  above,  and  h 
can  be  calculated.  If  the  value  of  h is  greater  than  the 
tolerance  value,  suitable  correction  for  the  assumed  K 
value  (i.e.,  K‘)  should  be  made,  and  the  whole  com- 
putation should  be  repeated.  For  stream  sanitation 
work,  a tolerance  value  of  0.01  for  h is  adequate. 

The  BODs  exerted  by  waste  samples  are  determined 
under  standardized  laboratory  conditions  at  20°C. 
Consequently,  the  values  of  ultimate  BOD,  K^,  and  K„ 
obtained  from  the  BOD  progression  data  pertain  to 
20°C.  To  assess  the  impacts  of  waste  discharges  at  tem- 
peratures other  than  20°C,  corrections  for  these  values 
must  be  applied.  The  following  equations  are  used  for 
temperature  correction  purposes. 


— K XI 

c(T)  '^c(20)  ^ ‘ ' 

[Eq  27] 

— K XI  007^“^®^ 

n(T)  *'n(20)  ^ i.W/ 

[Eq  28] 

Range  10  to  22°C 

n(T)  = K„ao)X  1.203X  0.877(^-22) 

[Eq29l 

Range  22  to  30°C 

'uc(T)  - Luc(20)  X (0.02T  + 0.6) 

[Eq  30] 

understood.  Various  theories  have  been  proposed  to 
explain  this  mechanism.  Some  of  these  theories  use 
parameters  like  surface  film  thickness,  surface  removal 
rate,  etc.;  however,  these  parameters  are  difficult  to 
measure  or  to  relate  to  common  hydraulic  variables. 
Therefore,  many  empirical  and  semi-empirical  equations 
have  been  proposed  to  predict  the  reaeration  rate.  The 
following  paragraphs  briefly  summarize  some  of  these 
conceptual  models  and  empirical  equations. 

Of  the  conceptual  models,  two  film  models  by  Lewis 
and  Whitman,  the  penetration  theory  model  by  Higbie, 
the  surface  renewal  theory  model  by  Dankwerts,  and 
the  film  penetration  model  by  Dobbins  are  well  known. 
All  these  models  contain  some  intangible  parameter 
which  cannot  be  directly  related  to  commonly  measured 
hydraulic  variables.  While  these  models  can  be  used  to 
explain  the  reaeration  process,  they  are  not  very  useful 
for  predicting  the  reaeration  rates  in  streams.  The 
following  prediction  equations  can  be  applied  more 
practically. 

Prediction  Equations 

O’Coiuior-Dobbins  Equation.  O’Connor  and  Dobbins 
defined  the  rate  of  reaeration,  Kj  to  the  base,  e,  in  terms 
of  the  rate  of  renewal  of  surface  film,  and  assumed 
that  the  best  estimate  was  given  by  the  ratio  of  the 
vertical  velocity  fluctuation  and  mixing  length.  In 
streams  in  which  there  is  pronounced  velocity  gradient, 
the  following  equation  for  nonisotropic  turbulence  was 
developed: 


where  T is  the  temperature  in  degrees  Celsius  of  the 
waste  stream  under  field  conditions,  is  the  value 
of  the  carbonaceous  deoxygenation  coefficient  at  T 
degrees  Celsius,  is  the  value  at  20°C,  and  so  on 
for  the  rest  of  the  parameters  defined  earlier. 


Kj  = 480D^®sO-^*H-i  25  [Eq  31] 

For  isotropic  turbulence  in  relatively  deep  channels 
where  there  is  no  pronounced  velocity  gradient. 


Reaeration  and  Photosynthesis 

Many  factors  govern  the  DO  content  of  rivers, 
streams,  and  other  water  bodies.  Oxygen  in  the  waters 
is  consumed  during  the  process  of  organic  materials 
stabilization,  oxidation  of  reduced  compounds  like 
ammonia,  ferrous  iron,  and  manganese  compounds, 
and  during  aquatic  plant  respiration.  Oxygen  is  replaced 
in  natural  water  bodies  largely  through  atmospheric 
reaeration  and  photosynthesis  of  aquatic  plants.  Evalua- 
tion of  the  river’s  DO  content  depends  on  the  rates  at 
which  the  processes  listed  above  occur. 

Even  though  atmospheric  reaeration  is  the  main 
source  of  oxygen  replenishment  in  natural  waters,  the 
mechanism  by  which  oxygen  is  absorbed  is  not  well 


Ki  = 1 27D^®sO-2Sh-i.s  [£q  32] 

where  D,„  is  the  coefficient  of  molecular  dispersion,  S 
is  the  channel  slope,  and  H is  the  mean  depth.  Isotrop- 
ic turbulence  was  assumed  arbitrarily  for  Chezy  co- 
efficients greater  than  17  and  nonisotropic  turbulence 
was  assumed  for  coefficients  less  than  17.  The  authors 
subsequently  concluded  that  the  latter  of  the  two 
equations  is  the  more  generally  applicable  and  recom- 
mended that  the  first  equation  be  ignored.^ 


*W.  E.  Dobbins,  “1964  BOD  and  Oxygen  Relationship  in 
Streams”  Journal  of  the  Sanitary  Engineering  Division,  Vol90, 
No.  SA3  (ASCE,  1964),  pp  53-78. 
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Churchiil-Elinoiv-Buckinghain  Equation.  The  empir- 
ical equation  proposed  by  Churchill,  et  al.,  is  based  on 
regression  analysis  of  the  extensive  field  data  from 
selected  stream  reaches  within  the  Tennessee  Valley 
Authority  rivers.  The  data  collected  by  these  investi- 
gators rank  as  one  of  the  best  efforts  in  gathering  field 
information  on  stream  reaeration  rates.  The  proposed 
empirical  equation  is 

Kj  = [Kq.^.H| 

in  which  Kj  is  the  reaeration  rate  coefficient  to  the  base, 
e,  day“';  V is  the  mean  velocity  of  flow  in  feet/sec; 
and  H is  the  mean  depth  of  flow  in  feet.  The  equation 
was  derived  with  a set  of  values  for  V and  H ranging, 
respectively,  from  1.85  to  5.00  ft/sec,  and  from  2.12 
to  11.41  ft.  Consequently,  use  of  this  equation  is 
restricted  to  the  values  of  mean  velocity  and  mean  depth 
that  fall  within  these  ranges. 

Isa'acs  and  Gaudy  Equation.  Isaacs  and  Gaudy  mea- 
sured reaeration  in  a circular  flume  in  which  the  inner 
and  outer  walls  could  be  driven  at  different  speeds  to 
simulate  stream  flow  conditions.  Multiple  regression 
analysis  on  their  data  yielded  the  equation 

Kj  = 7.05Vir^/2  [Eq  34] 

All  the  terms  are  as  defined  in  the  previous  equation. 

Langbien  and  Durum  Equation.  Langbien  and 
Durum  developed  an  empirical  equation  using  all  the 
published  data  on  atmospheric  reaeration.  The  pro- 
posed equation  is 

Kj  =7.63VH-'”  IEq35] 

All  the  terms  are  defined  in  the  Churchill-Elmore- 
Buckingham  equation. 

Several  other  semi-empirical  equations  predict  the 
reaeration  rate  coefficient  K} . Equations  proposed  by 
Krenkel  and  Orlob,  Thackston  and  Krenkel,  and  others, 
involve  determination  of  longitudinal  dispersion  in 
streams,  which  is  expensive,  both  in  terms  of  time  and 
manpower.  Prediction  equations  proposed  by  Church- 
hill,  ct  al.,  and  by  Langbien  and  Durum  arc  reliable  and 
acceptable  for  use.’ 


Temperature  Correction  for  Reaeration  Coefficient 
All  the  empirical  equations  given  above  are  for  deter- 
mining the  reaeration  coefficient  at  20°C.  However, 
when  the  stream  temperatures  are  other  than  20°C, 
corrections  must  be  applied.  A widely  accepted  cor- 
rection factor  for  the  reaeration  coefficient  is 

K2(t)  = Kj(jo)X  1.0241^-^®  [Eq36l 

where  and  are  the  values  of  reaeration 

coefficients  at  T°C  and  20°C,  respectively. 

The  reaeration  process  due  to  absorption  from  the 
atmosphere  is  considered  as  a first-order  reaction.  The 
rate  of  oxygen  absorption  is  dependent  on  the  DO 
deficit,  D in  mg/C  (which  is  defined  as  the  difference 
between  the  DO  saturation),  Cj  in  mg/C,  and  the  DO 
concentration  C in  mg/C  at  any  time,  t.  The  process  can 
be  expressed  as 


^=K,(C,-C)  = K,D  = -^ 

lEq  37] 

If  the  initial  DO  concentration  and  deficit  are  defined 
as  Cq  and  D^  at  time  t = 0,  the  solution  to  the  differ- 
ential Eq  37  is 

C,-C  = (C,-C„)e-'^7' 

[Eq  38] 

or 

D=  D„e-'^2' 

[Eq  39] 

Eqs  38  and  39  define  the  DO  deficit  with  time  when 
only  the  reaeration  process  is  operative  and  no  other 
process  like  deoxygenation, photosynthesis,  respiration, 
etc.,  is  affecting  the  system.  In  natural  water  bodies,  all 
these  processes  interact,  and  therefore  must  be  con- 
sidered as  coupled  systems. 

The  solubility  rate  of  oxygen  in  fresh  water  is  de- 
pendent primarily  on  the  water  temperature.  The 
saturation  concentrations  in  water  samples  were 
adequately  described'®  by  the  expression 

C,  = 14.652- 0.4 1022T 

+ 0.00799 IT’  - 0 000077774T'  (Eq  40) 
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’R.  E.  Rathben,  “Keaeration  Coefricients  of  Streams  - 
Stale  of  the  Art.”  J.  HvJr  Div.  Proc..  Vol  103,  No.  409 
(ASCE,  1977). 


'®T.  A.  Butts  and  V.  Kothandaraman,  “Fittinft  First  and 
Second  Order  BOD  Reaction  Equations  to  Stream  Data." 
Water  and  Sewage  Works.  Vol  1 1 7 (1970),  p 276. 
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where  Cj  is  the  saturation  concentration  in  nig/8  and  T 
is  the  water  temperature  in  degrees  Celsius.  Evaluation 
of  the  C,  term  at  any  specific  temperature  is  essential. 
Most  of  the  models  describing  the  DO/BOD  relation- 
ships evaluate  oxygen  deficits.  Thus,  knowing  the  DO 
saturation  concentrations  and  the  DO  deficits,  expected 
DO  concentrations  could  be  determined  by  finding  the 
difference  between  the  C,  value  and  the  DO  deficit 
value.- 

Oxygen  balance  studies  of  polluted  streams  usually 
produce  one  or  more  dissolved  oxygen  profiles.  Such 
a profile  is  referred  to  as  a dissolved  oxygen  sag  curve. 
The  analysis  of  a stream’s  capacity  to  assimilate  organic 
pollution  has  been  based  on  the  classical  theory  of 
Streeter  and  Phelps,*'  who  considered  only  two  mech- 
anisms affecting  DO  budget  in  streams;  (1)  the  bio- 
chemical oxidation  of  organic  matter,  and  (2)  the 
atmospheric  reaeration  of  river  water.  Since  their 
pioneering  efforts  several  modifications  have  been 
proposed. 

The  predictive  equations  developed  in  this  report 
for  BOD  and  DO  profiles  along  a stream  are  based  on 
the  following  assumptions: 

1 . Stream  flow  is  steady  and  uniform. 

2.  BOD  and  DO  are  uniformly  distributed  over 
each  cross  section,  so  that  the  equation  can  be  written 
in  one-dimensional  form. 

3.  The  process  for  any  river  reach  as  a whole  is  a 
steady-state  process,  with  the  conditions  at  every  cross 
section  remaining  unchanged  over  time. 


5.  The  rate  of  BOD  removal  is  proportional  to  the 
amount  of  BOD  present,  and  the  rate  of  reaeration  is 
proportional  to  the  oxygen  deficit. 

Factors  affecting  the  DO/BOD  relationship  in 
streams  considered  in  these  equations  are. 

1.  Oxygen  demand  by  carbonaceous  matter 

2.  Oxygen  demand  by  nitrogenous  matter 

3.  Oxygen  addition  due  to  atmospheric  reaeration 

4.  Oxygen  addition  due  to  photosynthesis. 

Other  factors  which  could  affect  the  oxygen  balance, 
such  as  biological  extraction  of  organic  matter,  benthal 
oxygen  demand,  and  deposition  or  scouring  of  organic 
matter  from  the  stream  beds,  are  ignored  as  having  a 
negligible  effect.  In  considering  the  oxygen  addition  to 
streams  because  of  photosynthesis.  Camp'*  used  a 
uniform  net  rate  of  oxygen  addition,  “a”  mg/?/day, 
despite  the  fact  that  there  is  no  oxygen  addition  from 
photosynthesis  during  periods  of  cloud  cover  or  dark- 
ness. However,  O’Connor'*  first  introduced  the  concept 
that  photosynthetic  oxygen  addition  could  be  treated 
as  a cyclic  function  with  zero  values  during  periods  of 
darkness  and  as  a specific  functional  value  during  day- 
light periods.  Both  concepts  are  discussed  in  the  follow- 
ing two  sections.  Although  the  O’Connor  and  Dobbins 
method  is  believed  to  be  more  accurate.  Camp’s  method 
requires  less  input  data.  It  is  left  to  the  user  to  evaluate 
the  importance  of  photosynthesis  and  to  determine 
which  method,  if  any,  should  be  used.  Figure  4 shows 
basic  elements  of  dissolved  oxygen  balance  models  in 
streams. 


4.  Bacterial  oxidation  of  carbonaceous  and  nitro-  Oxygen  Balance  Equations.  Case  I 
genous  matter  and  the  absorption  of  oxygen  from  the  The  rate  of  BOD  removal  in  a stream  on  the  basis  of 

atmosphere  are  first-order  reactions.  first-order  reaction  kinetics  is  given  by 


" H.  W.  Streeter  and  E.  B.  Phelps,  "A  Study  of  the  Pollution 
and  Natural  Purification  of  the  Ohio  Rivet,”  Public  Health 
Bulletin  146  (U.S.  Public  Health  Service,  1925). 

*^T.  R.  Camp,  Wafer and/ri/mpuri't/es (Reinhold  Publishing 
Corporation,  1963). 

'*W.  E.  Dobbins,  “1964  BOD  and  Oxygen  Relationship  in 
Streams.”  Journal  of  the  Sanitary  Engineering  Division,  Vol  90, 
No.  SA3  (ASCE,  1964),  pp  53-78. 

'^D.  J.  O’Connor  and  W.  E.  Dobbins  “Mechanisms  of 
Reaeration  in  Natural  Streams,”  rrarts.  ASCE,  Vol  123  (ASCE, 
1959). 


^ = -K,L,-K„L„  [Eq41| 

where 

L = total  BOD  remaining  at  any  time  t,  mg/S 
Lj  = carbonaceous  BOD  remaining  at  any  time  t, 
mg/e 


'*T.  R.  Camp. 
'*0’Connot  and  Dobbins. 
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Figure  4.  Basic  model  elcmenis. 


L„  = nitrogenous  BOD  remaining  at  any  time  t. 
mg/li 

Kj  = carbonaceous  deoxygenation  coefficient  to 
the  base  e,  per  day 

K„  = nitrogenous  deoxygenation  coefficient  to  the 
base  e.  per  day 

Integration  of  bq  41  with  appropriate  initial  conditions 
yields 


L = L e'"*'*^*  + L 


Ibq  42) 


where  is  the  ultimate  carbonaceous  demand  and 
is  tite  ultimate  nitrogenous  demand. 

If  there  is  a time  lag  (t^,)  for  the  onset  of  nitrogenous 
demand  from  the  carbonaceous  demand,  which  invari- 
ably is  the  case,  then  bq  42  takes  the  form 

L=L„,e-''c' + L„„e-'^n"-‘o)  [£^43] 

The  reaction  process  in  a stream  due  to  atmospheric 
absorption  is  also  a first-order  reaction  The  rate  process 
is  dependent  on  the  DO  deficit  (D  in  mg/C)  which  is 
defined  as  the  difference  between  the  DO  saturation 
concentration  (C,  in  mg/C)  and  the  observed  DO  con- 
centration (C  in  mg/C)  in  the  stream.  If  oxygen  addition 
from  photosyntliesis  is  considered  to  be  at  a uniform 
rate  of  "a"  mg/C/day.'^  tlien  the  rate  of  change  of  DO 
concentration  resulting  from  only  these  two  processes 
is  given  by 


’ T.  R.  Camp,  Miter  and  /is  /mpuritieti Reinhokl  Puhlishinf 
C'orptrratiun,  t963). 


^=Kj(C~C)  + a=KjD  + a = -^  (Eq44l 

Since  L,  the  BOD  remaining  to  be  satisfied,  is  measured 
in  terms  of  oxygen  consumed,  and  the  rate  of  oxygen 
consumption  equals  the  rate  at  which  the  BOD  is  satis- 
fied, then 


dt  dt 


IEq451 


Combining  the  rates  of  these  two  reactions  (Eqs  41  and 
44)  and  writing  the  equation  in  terms  of  DO  deficit 
yields  the  response  of  receiving  waters  to  a single  waste 
source.  The  differential  equation  is 


f + K„L„„e-'^"<'-o'-  KjD  - a 


(Eq  46| 


Solution  to  the  differential  Eq  46  is 


DeKjt  = + _121ILcK2t_Kn(l-to) 

K,-K,  K,-k/ 


+ constant 


[Eq471 


When  t = 0 (i.e.,  at  the  time  of  initial  observation, 
D = Dg  (initial  DO  deficit] ).  This  condition  yields 

Constant  = D„  - ^ -f  — 

IVj  IXj  Nj  Kj 


Substituting  this  value  uf  the  constant  in  Eq  47, 


rate  is  assumed  to  vary  with  the  amount  of  sunlight 
(reaching  a peak  at  noon  and  a zero  value  at  sunset  and 
sunrise),  then  this  source  of  oxygen  could  be  defined 
by  a periodic  function^  such  as  the  following: 


-j^(l-e-'^2')  + D„e-'^2'  (Eq481 

Eq  48  gives  a trace  of  the  DO  deficit  along  with  stream 
stretch.  Knowing  the  velocity  of  flow  and  the  distance 
of  the  point  of  interest  (i.e.,  the  point  at  which  DO 
concentration  must  be  determined)  from  the  initial 
point  of  reference,  the  elapsed  time  (t  in  days)  can  be 
evaluated.  The  DO  deficit  can  be  obtained  from  Eq  48, 
and  the  DO  concentration  at  the  desired  river  point  can 
then  be  deduced. 


P„,sin^t' 


when  0 < t'  < p 
when  1 — p < t' 


(Eq49| 


where  p is  the  period  of  oscillation  of  the  periodic 
function,  p^  is  the  amplitude  of  the  periodic  function, 
and  t is  the  time  of  flow  in  days  of  a mass  of  water 
since  sunrise  after  the  mass  of  water  has  entered  the 
head  end  of  the  river  reach  being  considered.  Figure  5 
illustrates  this  periodic  function. 


If  the  period  of  oscillation,  p,  is  taken  as  1 2 hours 
or  half  a day,  the  periodic  function  shown  in  Eq  49 
reduces  to 


Input  for  this  analytical  model  is  K^,  Kj,  K„, 
Lyn,  ty,  a,  and  D„.  Methods  for  evaluating  the  param- 
eters Kj,  K„,  Lu„,  Kj,  and  t„  have  been  discussed 
previously.  Oxygen  addition  due  to  photosynthesis  of 
aquatic  plants  (“a”)  must  be  determined  in  the  field 
using  the  light  and  dark  bottle  technique.'^  Most 
frequently,  the  photosynthesis  oxygen  addition  in 
flowing  waters  is  insignificant  and  can  be  ignored  in 
defining  the  oxygen  balance  equations.  It  is  included 
in  the  model  to  cover  those  cases  where  photosynthesis 
is  important. 

Qxygen  Balance  Equations,  Case  2 (O’Connor) 

The  drawback  of  the  Camp  approach  is  its  failure  to 
account  for  oxygen  addition  due  to  photosynthesis; 
realistically,  a uniform  rate  covering  the  daylight  and 
darkness  periods  is  considered,  even  though  no  oxygen- 
producing  photosynthesis  occurs  after  sunset.  This 
drawback  is  overcome  in  the  O’Connor  method.” 

Photosynthesis  depends  on  many  factors,  such  as  sun- 
light, temperature,  mass  of  algae  present,  and  nutrients. 
If  the  effects  of  these  factors  are  included  in  the  term 
P,,  which  represents  the  overall  rate  at  which  oxygen  is 
released  by  photosynthesis,  and  if  the  photosynthetic 


'*T.  Gaardner  and  H.  H.  Gran,  “Investigations  of  the  Pro- 
duction of  Plankton  in  the  Oslo  Fjord,”  Conseil  Permanent 
Internal.  Explor.  Mer.  Process.  Verhaux  Reunione,  Vol  42 
(1927),  pp  1-48. 

”d.  J.  O’Connor  and  W.  E.  Dobbins  “Mechanisms  of 
Reaeiation  in  Natural  Streams,”  Trans.  ASCE,  Vol  123  (ASCE, 
1959). 


P.'  = 


P^jSin  27rt' 
0 


when  0 < t'  < 1/2 

[Eq  50) 

when  l/2<t'<  1 


This  periodic  function  can  be  represented  quite 
accurately  by  the  first  three  terms  of  the  Fourier  Series, 
as  shown  below. 


P,’ — — -F -:^sin  27rt' — r-^cos4fft'  [EqSlJ 

‘ TT  2 JTT 

This  approach  of  expressing  the  diurnal  variations  in 
photosynthetic  activity  quantitatively  is  decidedly 
superior  to  averaging  them  over  the  entire  period,  wltich 
would  include  the  nighttime  flow  conditions. 

Using  this  expression  (Eq  51)  for  photosynthetic 
oxygen  addition,  the  oxygen  balance  equation  becomes 

+ K„L„„e-'^"<‘-‘o)  - K,D 

— P^(^-F  1/2  sin  27rt' — ^ cos  dirt') -F  R 

lEq52i 

where  R is  the  rate  of  respiration  (mg/8/day)  of  the 
aquatic  biota  and  all  the  other  terms  are  as  defined 
earlier. 


*0’Connor  and  Dobbins. 


e)  PERIOD  OF  OSCILLATION  P DAY 


TIME 


b)  PERIOD  OF  OSCILLATION  1/2  DAY 


SAM. 


I2N00N 


6RM. 

TIME— 


SAM. 


Figure  S.  The  periodic  function. 


The  solution  to  this  differential  equation  is 


[)e'<2'  = LJl!i_e(K2-K,.). 

K2  - K^, 

K,  - K„  rrK, 

P c''2' 

~ 2(Ki'+ 4;r^  ) 

2P  6*^2' 

+ — ^ :r  [4n  sin  rrt'  + K2  cos  4n’t'j 

37r  K?  + Ibrr^ 


Applying  the  initial  conditions  of  D = when  t = 0, 


*^c^uc 

Constant  = D — , 4-  — — 

•>  Kj  - K^.  Kj  - k/  ttKj 


^l+4V)''^'“"2'"~22rcos2rrr| 


3rr(KJ  + I6»r^) 


— (4irsin  rrr  + Kj  cos4ffr| 


lEq  54| 


R e'^2' 

+ + constant 

Kj 


[Eq531 


In  Eq  54,  the  term  r,  in  days,  is  the  elapsed  time  since 
sunrise  (6  AM)  for  the  initial  observation,  i.e.,  when 
t'  = 0,  t = T. 


■-.2  .^1  ■.  '.c-.  -.JL  . ♦.  . 


I 


I 


Substituting  tlie  value  uf  the  constant  in  Eq  S3,  the 
DO  deficit  equation  is  given  as: 


13= 


K,-K, 


K L 

+ .iL^"  . [e-K„(l-lo)_e-('^2*-Knto)] 


K2-K„ 


^ , R V , 

-;^(l-e-'‘2')  + i^(l-e-'^i') 


2(K|  + 4jr’) 


[(K2  sin2irt' — 27r  cos2fft') 


— e sin2ffr-2ir  cos27rT)] 

2P„ 


— e”  *^2' (47r  sin7rT+K2Cos47rT)| 

+ D„e-'‘J'  [Eq  55) 


Once  the  DO  deficit  is  assessed  at  any  point  in  a given 
river  reach,  the  dissolved  oxygen  can  be  estimated  as 
discussed  earlier. 


If  the  BOD  resulting  from  the  oxidation  of  nitro- 
genous matter,  oxygen  addition  from  photosynthesis, 
and  aquatic  biota  respiration  are  considered  insignificant 
compared  to  the  two  important  processes  affecting  the 
oxygen  balance  (the  BOD  due  to  the  oxidation  of  car- 
bonaceous matter  and  the  atmospheric  reaeration), 
then  the  DO  deficit  equations  (Eqs  48  and  55)  reduce 
to  the  classical  form  given  by  Streeter  and  Phelps. 


D = n'^4^(e-'^c'_e-K2<)+D„e-'^j'  (Eq  56) 

IS.2  IS-c 


H.  W Streeter  and  E.  B.  Phelps,  "A  Study  of  the  Pollution 
and  Natural  Purification  of  the  Ohio  River."  Public  Health 
Bulletin  l46iV.S.  Public  Health  Service,  I92S). 


3 USE  OF  THE  DOSAG  MODEL 


Batch  Mode 

DOSAG  combines  the  three  analysis  functions  into 
one  computer  program.  Ten  output  options  are  avail- 
able to  the  user,  as  described  in  Table  1.  To  obtain  out- 
put, the  user  should  first  complete  the  input  form 
(Figure  6),  supplying  all  the  required  input  data.  Table 
2 defines  the  required  input  parameters,  and  Figure  7 
illustrates  the  format  for  submitting  the  paired  observa- 
tions data  required  for  analysis  options  2 and  3.  The 
input  form  should  then  be  mailed  to: 


US  Army  CERL/EN 
ATTN:  Bob  Riggins 
P.O.  Box  4005 
Champaign,  Illinois  61820 


CERL  will  then  put  the  data  into  the  format  neces- 
sary for  running  it  througli  the  computer  and  mail  the 
resulting  output  to  the  user. 


Figures  8 througli  1 1 are  examples  of  the  types  of 
output  available  from  DOSAG.  The  test  data  are  fic- 
tional and  were  developed  to  illustrate  as  many  options 
as  possible. 


ConverMtional  Mode 

A computer  program  operating  in  conversational 
mode  prompts  the  user  for  input  data  and  option  selec- 
tion. The  user  is  given  a list  of  nine  options  from  which 
to  choose;  he/she  types  “1”  if  the  option  is  chosen,  or 
“0”  if  it  is  not.  Following  are  the  available  options. 


OPTION  1-USE  OBSERVED  DATA  TO  CALC 
COEFF  FOR  SINUSOIDAL  ANNUAL  WATER  TEMP 
CURVE  IF  DESIRED  TYPE  1 IF  NOT  0 


OPTION  2-USE  OBSERVED  DATA  TO  CALC 
CARBONACEOUS  DEOXYGENATION  COEFF,  VIA 
BOD,  TYPE  1 OR  0 


OPTION  3-USE  OBSERVED  DATA  TO  CALC 
NITROGENOUS  DEOXYGENATION  COEFF.  VIA 
BOD,  TYPE  1 OR  0 


OPTION  4-CALC  DISSOL  OXYGEN  SAG  FOR  A 
SINGLE  REACH  OF  STREAM  WITH  ONE  WASTE 
DISCHARGE.  VIA  DO,  TYPE  1 ORO 
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Table  1 

DOSAG  Options 


Table  2 

Parameter  Definitions 


Option  Description  Parameter 


1 

Calculate  a dissolved  oxygen  profile  for  a single 

1.1 

UPSTO 

reach  of  stream  with  one  waste  inflow,  considering 

1.2 

UPSTDO 

only  carbonaceous  deoxygenation  demand. 

1.3 

UPSBODC 

2 

Include  nitrogenous  deoxygenation  demand  in  the 

1.4 

EFFQ 

calculations  for  dissolved  oxygen  profile. 

1.5 

EFFDO 

1.6 

EFFBODC 

3 

Use  Churchill  method  for  reacration.  If  this  option 
is  not  specified,  the  Langbien-Durum  equation  will 

1.7 

WIDTH 

be  used. 

1.8 

DEPTH 

1.9 

DISTC 

4 

Include  photosynthesis  deoxygenation  in  the  cal- 
culations for  dissolved  oxygen  profile  using  Camp’s 

1.10 

KSUBC 

method. 

1.11 

TEMP 

5 

Include  photosynthesis  deoxygenation  using 

2.1 

UPSBODN 

O’Connor’s  method. 

2.2 

EFFBODN 

6 

Use  sinusoidal  curve  coefficients  (which  have  been 
previously  determined)  to  calculate  water  tempera- 

2.3 

KSUBN 

ture  for  use  in  determining  dissolved  oxygen  profile. 

2.4 

TZERO 

7 

Use  with  option  10  to  specify  a day  for  which  a 
temperature  value  is  to  be  calculated.  This  temper- 
ature value  will  be  used  if  option  1 is  selected. 

4.1 

PSUBM 

8 

Use  observed  data  to  calculate  a carbonaceous 

4.2 

RESP 

deoxygenation  rate  coefficient  and  ultimate  bio- 
chemical oxygen  demand.  (If  option  1 is  also 
chosen,  these  parameters  will  become  input  data 

4.3 

TALL 

to  the  calculations  for  dissolved  oxygen  profile.) 

5.1 

A 

9 

Use  observed  data  to  calculate  a nitrogenous  deoxy- 
genation  rate  coefficient  and  ultimate  biochemical 
oxygen  demand.  (If  options  1 and  2 are  also 
chosen,  these  parameters  will  become  input  data  to 

6.1 

P 

the  calculations  for  dissolved  oxygen  profile.) 

6.2 

0 

10 

Use  observed  data  to  calculate  coefficients  for 

6.3 

TBAR 

sinusoidal  annual  water  temperature  curve.  (If 
option  1 is  also  chosen,  a value  for  water  tempera- 

6.4 

DAY 

ture  will  become  input  data  for  calculations  of 
dissolved  oxygen  profile.) 

7.1 

DAY 

8.1 

EKSUBC 

Paired  Observations 


9.1  EKSUBN 


Paired  Observations 
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Definition 

Upstream  flow  (cfs) 

Upstream  dissolved  oxygen  (mg/S) 
Upstream  carbonaceous  deoxygena- 
tion demand  (mg/e) 

Effluent  flow  (cfs) 

Effluent  dissolved  oxygen  (mg/e) 
Effluent  carbonaceousdeoxygenation 
demand  (mg/e) 

Average  stream  width  (ft) 

Average  stream  depth  (ft) 

Length  of  the  stream  reach  (miles) 
Carbonaceous  deoxygenation  rate 
coefficient 

Water  temperature  (°C) 

Upstream  nitrogenous  deoxygenation 
demand  (mg/e) 

Effluent  nitrogenous  deoxygenation 
demand  (mg/e) 

Nitrogenous  deoxygenation  rate 
coefficient 

Time  lag  for  nitrogenous  deoxygena- 
tion demand  (days) 

Amplitude  of  the  photosynthesis  peri- 
odic function 

Rate  of  respiration  (mg/e/day) 

Time  of  flow  (days)  since  sutuise  after 
the  mass  of  water  has  entered  the 
head  end  of  a river  reach 

Oxygenation  rate  due  to  photosyn- 
thesis (mg/fi/day) 

Sinusoidal  water  temperature  curve 
coefficient 

Sinusoidal  water  temperature  curve 
coefficient 

Average  water  temperature  (°C) 

Day  of  year  (Julian) 

Day  of  year  (Julian) 

Estimated  carbonaceous  deoxygena- 
tion rate  coefficient 

Observed  values  of  carbonaceous 
BOD  (mg/S)  over  time  (days)  deter- 
mined by  laboratory  analysis 

Estimated  nitrogenous  deoxygenation 
rate  coefficient 

Observed  values  of  nitrogenous  BOD 
(mg/t)  over  time  (days)  determined 
by  laboratory  analysis 


1 


,1 


DIRECTIONS  - Check  the  options  required  DATE  

and  provide  the  necessary  input  NAME 

information.  Mail  the  completed  ADDRESS 

input  form  to:  US  Army  CERL/EN  

P.  0.  Box  4005  PHONE 


Champaign,  IL  61820 


Notes:  1.  Options  2-7  require  option  1 . 

2.  Option  7 also  requires  option  10. 

3.  Values  1.6  and  1.10  can  be  furnished  by  option  8. 

4.  Value  1.11  can  be  furnished  by  option  6 or  7. 

5.  Values  2.2  and  2.3  can  be  furnished  by  option  9. 

6.  Paired  observations  data  should  be  submitted  on  a 
separate  data  sheet  such  as  that  illustrated  in 
Figure  7 of  the  user  manual. 


1.1  UPSTO 

1.4  EFFO 

1.7  WIDTH 

1.2  UPSTDO 

1.5  EFFDO 

1.8  DEPTH 

1.3  UPSBODC 

1 .6  EFFBODC 

1.9  DISTC 

1.10  KSUBC 

1.11  TEMP 

Figure  6.  Input  form  for  DOSAG. 


□ Option  8 


□ Option  9 


□ Option  10 


9.1  EKSUBN 

Paired  Observations: 

9.2  times 

9.3  BODN  values 


□ Option  2 


□ Option  4 


□ Option  6 


2.1  UPSBODN 

2.2  EFFBODN 

2.3  KSUBN 

2.4TZERO 


4.1  PSUBM 

4.2  RESP  _ 

4.3  TAU 


6.1  P 

6.2Q 

6.3  TBAR. 

6.4  DAY 


□ Option  3 


□ Option  5 


□ Option  7 


□ Option  1 


8.1  EKSUBC 

Paired  Observations: 

8.2  times 

8.3  BODC  values 


Paired  Observations 

10.1  Days 

10.2  Temp,  values 


21 


ENVIRONMENTAL  DIVISION 

CONSTRUCTION  ENGINEERING  RESEARCH  LABORATORY 
UNITED  STATES  ARMY  CORPS  OF  ENGINEERS 


! 


ANALYTICAL  TECHNIQUES  FOR  ENVIRONMENTAL  IMPACT  ANALYSIS 
WATER  QUALITY  ANALYSIS  SYSTEMS 
DOSAG  COMPONENT 
SEPTEMBER  1977  VERSION 


OUTPUT  FOR: 

BAYOU  ANACOCO  - FORT  POLK,  LA 


ANALYSIS  OPTIONS  SELECTED; 

T — USE  OBSERVED  DATA  TO  CALCULATE  COEFFICIENTS 
^ FOR  SINUSOIDAL  ANNUAL  WATER  TEMPERATURE  CURVE  AND,  IF  DO 

I IS  TO  BE  COMPUTED  CALCULATE  WATER  TEMPERATURE  FOR  A 

GIVEN  DAY  OF  THE  YEAR  VIA  STEMPl 

C — USE  OBSERVED  DATA  TO  CALCULATE  CARBONACEOUS 
DEOXYGENATION  COEFFICIENT  VIA  BOD 

N — USE  OBSERVED  DATA  TO  CALCULATE  NITROGENOUS 
DEOXYGENATION  COEFFICIENT  VIA  BOD 

D — CALCULATE  DISSOLVED  OXYGEN  SAG  FOR  A SINGLE 
REACH  OF  STREAM  WITH  ONE  WASTE  DISCHARGE  VIA  DO.  IF  C 
? IS  SPECIFIED,  USE  COMPUTED  KSUBC.  IF  Y IS  SPECIFIED 

CONSIDER  NITROGENOUS  DEOXYGENATION;  AND,  IF  N IS 
SPECIFIED,  USE  COMPUTED  KSUBN. 

A — USE  CHURCHILL'S  METHOD  FOR  REAERATION: 
OTHERWISE,  USE  LANGBERN-DURUM  EQUATION. 

: P — INCLUDE  PHOTOSYNTHESIS  USING  CAMP'S  METHOD 

Y — INCLUDE  NITROGENOUS  DEOXYGENATION  IN  THE  DO 
SAG  CALCULATION. 


Figure  8.  Example  DOSAG  output-options  selection. 


■i 
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SINUSOIDAL  CURVE  FITTING  FOR  OBSERVED  WATER  TEMPERATURES  BY  LEAST  SQUARE  FIT 


DAYS  FROM  INITIAL  PERIOD 

257.0 

165.0 

163.0 

182.0 
226.0 

253.0 

324.0 

13.0 

43.0 

78.0 

113.0 

148.0 

336.0 

295.0 

314.0 

363.0 

19.0 

46.0 

74.0 

116.0 

137.0 

165.0 

228.0 

207.0 

196.0 

269.0 

196.0 


OBSERVED  TEMPERATURE 

26.0 

27.0 

26.0 

27.0 

28.0 

24.0 

13.0 

7.0 

11.0 

16.0 

24.0 

26.0 

20.0 

22.0 

15.0 

18.0 

10.0 

13.0 

20.0 

19.0 

24.0 

28.0 

26.0 

23.0 

24.0 

23.0 

23.0 


AMPLITUDE  =4.2 
STAND.  ERROR  =3.7 
STAND.  DEV.  = 6.0 
COEFF.  of  CORR.  = .794 

CURVE  COEFF.  AND  AVE.  ANNUAL  TEMP 
A = -1.934 
Q = -3.698 
TBAR  = 20.852 

CALCULATED  TEMP  FOR  DAY  95.  IS  19.2 

Figure  9.  Example  DOSAG  output-temperature. 


COMPUTED  TEMPERATURE 

23.8 

23.8 

23.7 
24.5 

24.9 

24.0 

19.3 

16.8 

16.8 

18.1 

20.4 

22.8 

18.5 
21.4 

20.0 

17.2 

16.7 

16.9 

17.9 

20.6 

22.1 

23.8 

24.8 

25.0 

24.9 

23.1 

24.9 


CARBONACEOUS  BOD  PROGRESSION  INPUT  DATA 


NITROGENOUS  BOD  PROGRESSION  INPUT  DATA 


COMPUTED  TOTAL  BOD  PROGRESSION 


TIME  IN  DAYS 


BOD  SATISFIED 


.990 

2.420 

1.860 

4.190 

2.950 

5.600 

4.850 

7.020 

6.820 

8.250 

8.860 

9.230 

CARBONACEOUS  RATE 

COEFF.  = .279 

CARB.  ULT.  BOD  = 

9.834 

TIME  IN  DAYS 

BOD  SATISFIED 

3.000 

2.420 

4.000 

4.190 

5.000 

5.600 

7.000 

7.020 

NITROGENOUS  RATE 

COEFF.  = .333 

NITR.  ULT.  BOD  = 

8.703 

TIME  IN  DAYS 

TOTAL  BOD  SATISFIED 

0.000 

0.000 

.500 

1.282 

1.000 

2.398 

1.500 

3.367 

2.000 

4.210 

2.500 

4.944 

3.000 

5.581 

3.500 

6.136 

4.000 

6.618 

4.500 

7.037 

5.000 

7.402 

5.500 

7.719 

6.000 

7.995 

6.500 

8.235 

7.000 

8.443 

7.500 

8.625 

8.000 

8.782 

8.500 

8.919 

9.000 

9.039 

9.500 

9.142 

10.000 

9.232 

Figure  10.  Example  DOSAG  output-rate  coefficients 
and  ultimate  biochemical  oxygen  demand. 
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CRITERIA  FOR  COMPUTATION  OF  DO  DEFICIT 


UPSTREAM 

FLOW= 

21.000 

EFFLUENT  FLOW= 

4.800 

UPSTREAM 

BODC= 

2.000 

EFFLUENT  BODC= 

9.834 

UPSTREAM 

D0= 

7.000 

EFFLUENT  D0= 

.700 

WIDTH= 

35.000 

DEPTH= 

3.000 

KSUBC= 

.269 

TEMP= 

19.155 

DISTANC= 

50.000 

UPSTREAM 

BODN=  2.000  EFFLUENT  BODN= 

8.703 

KSUBN= 

.308 

TZERO= 

2.000 

REAREATION  RATE  COEFFICIENT=  .435 

A-VALUE  FOR  CAMP'S  PHOTOSYNTHESIS  IS,  .100 


DISSOLVED  OXYGEN  PROFILE  FOR  DISTANCE  OF  50  00  MILES 


IN  DAYS 

DISTANCE  IN  MILES 

DO.  DEFICIT 

DISS.  OXYGEN 

0.00 

0.00 

3.35 

5.83 

.20 

.80 

3.13 

6.05 

.40 

1.61 

2.93 

6.25 

.60 

2.41 

2.74 

6.44 

.80 

3.22 

2.56 

6.62 

1.00 

4.02 

2.39 

6.79 

1.20 

4.82 

2.23 

6.95 

1.40 

5.63 

2.09 

7.09 

1.60 

6.43 

1.95 

7.23 

1.80 

7.24 

1.82 

7.36 

2.00 

8.04 

1.70 

7.48 

2.20 

8.85 

3.39 

5.79 

2.40 

9.65 

3.31 

5.87 

2.60 

10.45 

3.22 

5.96 

2.80 

11.26 

3.13 

6.05 

3.00 

12.06 

3.03 

6.15 

3.20 

12.87 

2.93 

6.25 

3.40 

13.67 

2.83 

6.35 

3.60 

14.47 

2.73 

6.45 

Figure  1 1.  lixaiuple  DOSAG  output-DO  deficit.  ^ 
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TIME  IN  DAYS 


DISTANCE  IN  MILES 


DO.  DEFICIT 


DISS.  OXYGEN 


I- 


3.80 

15.28 

2.62 

6.55 

4.00 

16.08 

2.52 

6.66 

4.20 

16.89 

2.42 

6.76 

4.40 

17.69 

2.32 

6.86 

4.60 

18.50 

2.22 

6.96 

4.80 

19.30 

2.12 

7.06 

5.00 

20.10 

2.03 

7.15 

5.20 

20.91 

1.93 

7.24 

5.40 

21.71 

1.84 

7.34 

5.60 

22.52 

1.75 

7.42 

5.80 

23.32 

1.67 

7.51 

6.00 

24.12 

1.58 

7.59 

6.20 

24.93 

1.50 

7.68 

6.40 

25.73 

1.43 

7.75 

6.60 

26.54 

1.35 

7.83 

6.80 

27.34 

1.28 

7.90 

7.00 

28.15 

1.21 

7.97 

7.20 

28.95 

1.14 

8.04 

7.40 

29.75 

1.07 

8.11 

7.60 

30.56 

1.01 

8.17 

7.80 

31.36 

.95 

8.23 

8.00 

32.17 

.89 

8.29 

8.20 

32.97 

.84 

8.34 

8.40 

33.77 

.79 

8.39 

8.60 

34.58 

.74 

8.44 

8.80 

35.38 

.69 

8.49 

9.00 

36.19 

.64 

8.54 

9.20 

36.99 

.60 

8.58 

9.40 

37.80 

.56 

8.62 

9.60 

38.60 

.52 

8.66 

9.80 

39.40 

.48 

8.70 

10.00 

40.21 

.44 

8.74 

10.20 

41.01 

.41 

8.77 

10.40 

41.82 

.37 

8.81 

10.60 

42.62 

.34 

8.84 

10.80 

43.42 

.31 

8.87 

11.00 

44.23 

.28 

8.90 

11.20 

45.03 

.26 

8.92 

11.40 

45.84 

.23 

8.95 

11.60 

46.64 

.21 

8.97 

11.80 

47.45 

.18 

9.00 

12.00 

48.25 

.16 

9.02 

12.20 

49.05 

.14 

9.04 

12.40 

49.86 

.12 

9.06 

12.44 

50.00 

.12 

9.06 

Figure  II.  (con’t) 
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OPTION  5-USE  CHURCHILL’S  METHOD  FOR 
REAERATION  OTHERWISE  USE  LANGBIEN- 
DURUM,  TYPE  1 OR  0 

OPTION  6-INCLUDE  PHOTOSYNTHESIS  USING 
CAMPS  METHOD,  TYPE  1 IF  DESIRED,  IF  NOT 
TYPEO 

OPTION  7- INCLUDE  PHOTOSYN  USING  O’CON- 
NOR’S METHOD  TYPE  1 OR  0 

OPTION  9-COMPUTE  TEMP  GIVEN  THE  EQUA- 
TION COEFF,  MEAN  ANN  TEMP,  + DAY  OF 
YEAR,  TYPE  I ORO 

The  user  should  first  obtain  the  required  input  data 
for  the  options  he/she  will  select. 

The  model  can  then  be  accessed  through  telephone 
hook-up,  using  a remote  terminal. 

The  interactive  user  may  access  the  system  by  com- 
pleting four  tasks: 

(1)  Setting  remote  terminal  settings.  (Full  duplex, 
30  cps,  online) 

(2)  Dialing  the  appropriate  computer  number.  Two 
numbers  are  presently  available:  (Commercial:  217/ 
333-1587  and  FTS:  957-1587) 

(3)  Logging  into  the  computer  by  typing  the  user 
name  and  the  appropriate  password.  Following  is  an 
example  of  user  response. 

CAC  Network  Unix  System 
log  in:  riggins 

Password : 

In  response  to  the  “password”  request,  the  user 
supplies  the  appropriate  response.  (The  password  will 
not  be  printed  on  the  output  page.)  The  user  may 
obtain  a password  by  contacting  CERL/EN;  ATTN: 
Bob  Riggins,  P.O.  Box  4005,  Champaign,  IL  61820. 

System  messages  will  appear  next,  as  in  the  follow- 
ing example. 

15  Apr  76  UNIX  rpO  dump  will  be  done  between 

1700-1800  tonight. 


The  System  will  go  down  at  1655;  please  be  off 
by  then. 

New  System  testing  periods. . .run  bullboard. 

Issty:  not  found 

No  mail 

This  account  has  been  changed  to  accommodate 

added  systems. 

Type  “dosag”  to  get  into  CERL’s  programs. 

These  messages  are  only  for  relaying  system  infor- 
mation to  users.  Any  changes  in  CERL’s  ETIS  approach 
will  be  identified  here.  The  user  should  wait  for  a 
response  from  the  computer  before  entering  additional 
commands.  This  signal  will  be  the  character 

(4)  Logging  into  the  DOSAG  model: 

% dosag 

The  following  are  program  prompts  as  they  appear 
in  the  options  in  the  DOSAG  model  program. 

OPTION  /-USE  OBSERVED  DATA  TO  CALC 
COEFF  FOR  SINUSOIDAL  ANNUAL  WATER  TEMP 
CURVE  IF  DESIRED  TYPE  1 IF  NOT  0 

Option  1 calculates  the  coefficient  for  sinusoidal 
annual  water  temperature  curve.  If  this  option  is  desired, 
the  correct  user  response  is  1 . Option  1 requires  that 
the  user  provide  responses  to  three  data  requests. 

1.  Enter  the  number  of  observations  for  the  year 
being  analyzed. 

2.  Enter  the  Julian  number  for  each  day  for  which  a 
temperature  observation  is  used.  When  finished,  type 
999. 

3.  Enter  the  observed  temperature  corresponding  to 
each  given  Julian  number.  When  finished,  type  999. 

Data  requests  2 and  3 may  require  more  than  one 
response.  If  so,  the  user  should  enter  a different  value 
on  each  line,  depressing  the  carriage  return  after  each 
value,  and  typing  999  when  all  values  are  entered.  The 
user  should  insure  that  the  sequence  of  temperature 
entries  corresponds  to  the  sequence  of  day  values. 

OPTION  2-USE  OBSERVED  DATA  TO  CALC 
CARBONACEOUS  DEOXYGENATION  COEFF, 

VIA  BOD,  TYPE  1 OR  0 


Option  2 calculates  the  carbonaceous  deuxygenation 
coefficient.  If  this  option  is  desired,  the  correct  user 
response  is  1.  Option  2 requires  the  user  to  provide 
responses  to  five  data  requests. 

1 . Enter  the  total  number  of  BOD  observations  for 
carbonaceous  demand. 

2.  Enter  initial  estimate  for  carbonaceous  rate  co- 
efficient. 

3.  Enter  TZERO. 

4.  Enter  each  day  of  carbonaceous  demand  observa- 
tion in  numerical  sequence.  When  finished,  type  999. 

5.  Enter  each  BOD  carbonaceous  demand  observa- 
tion in  sequence  corresponding  to  day  of  observation. 
When  finished,  type  999. 

Data  requests  4 and  5 may  require  more  than  one 
response.  If  so,  the  user  enters  a different  value  on  each 
line,  depressing  the  carriage  return  after  each  value,  and 
typing  999  when  finished. 

OPTION  J-USE  OBSERVED  DATA  TO  CALC 
NITROGENOUS  DEOXYGENATION  COEFF, 

VIA  BOD,  TYPE  1 OR  0 

Option  3 calculates  the  nitrogenous  deoxygenation 
coefficient  and  may  not  be  selected  without  first  select- 
ing option  2.  Option  3 requires  that  the  user  provide 
responses  to  nine  data  requests. 

1.  Enter  total  number  of  BOD  observations  for 
carbonaceous  demand. 

2.  Enter  total  number  of  BOD  observations  for 
nitrogenous  demand. 

3.  Enter  carbonaceous  rate  coefficient. 

4.  Enter  nitrogenous  rate  coefficient. 

5.  Enter  TZERO. 

6.  Enter  day  of  carbonaceous  demand  observation. 
Wlien  finished,  type  999. 

7.  Enter  BOD  carbonaceous  demand.  When  finished, 
type  999. 


8.  Enter  day  of  nitrogenous  demand  observation. 
When  finished,  type  999. 

9.  Enter  nitrogenous  demand.  When  finished,  type 
999. 

Data  requests  6,  7,  8,  and  9 may  require  more  than 
one  response.  If  so,  the  user  enters  a different  value  on 
each  line,  depressing  the  carriage  return  after  each 
value,  and  typing  999  when  finished. 

OPTION  4-Q.k\£  DISSOL OXYGEN  SAG  FOR  A 
SINGLE  REACH  OF  STREAM  WITH  ONE  WASTE 
DISCHARGE,  VIA  DO,  TYPE  I OR  0 

Option  4 calculates  the  dissolved  oxygen  sag  for  a 
single  reach  of  stream  with  one  waste  discharge.  This 
option  requires  that  the  user  provide  responses  to  1 1 
data  requests. 

1.  Enter  effluent  BODC. 

2.  Enter  KSUBC. 

3.  Enter  upstream  flow. 

4.  Enter  effluent  flow. 

5.  Enter  upstream  BODC. 

6.  Enter  upstream  DO. 

7.  Enter  effluent  DO. 

8.  Enter  stream  width  (in  feet). 

9.  Enter  stream  depth  (in  feet). 

10.  Enter  stream  temperature  (in  °C). 

1 1 . Enter  distance  (in  miles). 

OPTION 5-\iSE  CHURCHILL’S  METHOD  FOR 
REAERATION  OTHERWISE  USE  LANGBIEN- 
DURUM,TYPE  1 ORO 

A “1”  response  to  Option  5 enables  the  use  of 
Churchill’s  method  for  reaeration;  a “0”  response 
provides  the  Langbien-Durum  method.  This  option 
is  used  in  conjunction  with  option  4 and  may  not  be 
used  unless  a “1”  response  has  been  given  to  option  4. 
There  are  no  additional  data  requests  for  option  4. 
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OPTION  6- INCLUDE  PHOTOSYNTHESIS 
USING  CAMP’S  METHOD,  TYPE  I 
IF  DESIRED,  IF  NOT  TYPE  0 

Option  6 output  includes  the  photosynthesis  cal- 
culations using  Camp’s  method  for  finding  dissolved 
oxygen  sag.  This  option  may  only  be  used  in  conjunc- 
tion with  option  4 and  may  not  be  used  unless  a “1” 
response  has  been  given  to  option  4.  In  addition  to 
the  data  requests  for  option  4,  the  following  input 
is  also  required. 

Enter  A-Value  for  Camp  Photosynthesis 

OPTION  7- INCLUDE  PHOTOSYN  USING 
O’CONNOR’S  METHOD,  TYPE  1 OR  0 

Option  7 output  includes  the  photosynthesis  calcu- 
lations using  O’Connor’s  method  for  finding  dissolved 
oxygen  sag.  This  option  may  only  be  used  in  conjunc- 
tion with  a “1”  response  to  option  4.  In  addition  to 
the  data  requests  for  option  4,  the  following  input  is 
also  required; 

1.  Enter  amplitude. 

2.  Enter  rate  of  respiration. 

3.  Enter  TAU. 

O/T/OA  5- INCLUDE  NITROG  DEOXYGENATION 
IN  DOSAG  CALC,  TYPE  1 OR  0 

Option  8 output  includes  nitrogenous  deoxygenation 
calculations  for  dissolved  oxygen  sag.  This  option  may 
only  be  used  in  conjunction  with  a “1”  response  to 
option  4.  In  addition  to  the  data  requests  for  option  4, 
the  following  input  is  also  required; 

1.  Enter  upstream  BODN. 

2.  Enter  effluent  BODN. 

3.  Enter  KSUBN. 

4.  Enter  TZERO. 

OPTION  9-COMPUTE  TEMP  GIVEN 
THE  EQUATION  COEFF,  MEAN  ANN 
TEMP,  + DAY  OF  YEAR,  TYPE  1 OR  0 

Option  9 computes  the  temperature  for  a specific 
day  when  the  user  provides  the  equation  coefficients. 


mean  annual  temperature,  and  day  of  year.  This  option 
may  not  be  selected  separately,  but  may  be  used  in 
conjunction  with  a “ 1 ” response  to  either  options  1 or 
2.  In  addition  to  the  data  requests  for  options  1 or  2, 
the  following  input  is  also  required; 

1.  Enter  curve  coefficient  P. 

2.  Enter  curve  coefficient  Q. 

3.  Enter  average  annual  temperature  (in  °C). 

4.  Enter  the  Julian  day  for  which  temperature  will 
be  calculated. 

The  examples  beginning  on  page  31  show  the 
procedures  for  using  options  4 and  2,  respectively. 
Output  for  each  example  is  shown. 


4 USE  OF  OUTPUT 


Output  from  the  model  can  be  used  to  quantify  the 
impact  on  a stream’s  dissolved  oxygen  resources  result- 
ing from  an  Army  activity  which  introduces  to  the 
stream  material  exerting  a biological  oxygen  demand. 
Options  within  the  model  provide  intermediate  calcula- 
tions of  values  for  parameters  such  as  water  temperature, 
deoxygenation  rate  coefficients,  and  ultimate  biological 
oxygen  demand  which  may  be  needed  to  compute  the 
DO  sag  curve.  The  model  is  applicable  fora  single-reach, 
single-point-source  pollutant  inflow  situation.  The  flow 
conditions  must  correspond  to  a steady-flow-state. 

The  resulting  DO  sag  curve  profile  can  be  used  to 
determine  the  following; 

1 . The  minimum  DO  level 

2.  The  stream  location  where  DO  is  a minimum 

3.  The  stream  location  where  DO  level  recovers  to 
specified  levels  (saturation  level  or  quality  target  level). 

Analysis  of  this  information  will  indicate  the  magni- 
tude of  impact  on  the  stream’s  DO  resources. 
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Example  1 


% DOSAG 

OPTION  1— USE  OBSERVED  DATA  TO  CALC  COEFF 
FOR  SINUSOIDAL  ANNUAL  WATER  TEMP  CURVE 
IF  DESIRED  TYPE  1 IF  NOT  0 
0 

OPTION  2— USE  OBSERVED  DATA  TO  CALC  CARBONACEOUS 
DEOXYGENATION  COEFF.  VIA  BOD.  TYPE  1 OR  0 
0 

OPTION  3-USE  OBSERVED  DATA  TO  CALC  NITROGENOUS 
DEOXYGENATION  COEFF.  VIA  BOD.  TYPE  1 OR  0 
0 

OPTION  4-CALC  DISSOL  OXYGEN  SAG  FOR  A SINGLE 
REACH  OF  STREAM  WITH  ONE  WASTE  DISCHARGE  VIA  DO.  TYPE  1 OR  0 
1 

OPTION  5--USE  CHURCHILL'S  METHOD  FOR  REAERATION 
OTHERWISE  USE  LANGBERN-DURUM.  TYPE  1 OR  0 
0 

OPTION  6— INCLUDE  PHOTOSYNTHESIS  USING  CAMPS  METHOD. 

TYPE  1 IF  DESIRED,  IF  NOT  TYPE  0 
0 

OPTION  7-INCLUDE  PHOTOSYN  USING  O'CONNOR'S  METHOD 
TYPE  1 OR  0 
1 

OPTION  8-INCLUDE  NITROG  DEOXYGENATION  IN  DOSAG  CALC 
TYPE  1 OR  0 
0 

OPTION  9— COMPUTE  TEMP  GIVEN  THE  EQUATION  COEFF,  MEAN 
ANN  TEMP,  + DAY  OF  YEAR.  TYPE  1 OR  0 
0 

ENTER  EFFLUENT  BODC 
2 

ENTER  KSUBC 
2 

ENTER  UPSTREAM  FLOW 
2 

ENTER  UPSTREAM  BODC 
3 

ENTER  UPSTREAM  DO 


ENTER  EFFLUENT  DO 

3 

ENTER  STREAM  WIDTH 
5 

ENTER  STREAM  DEPTH 
2 

ENTER  STREAM  TEMP(IN  C) 

21 

ENTER  DISTANCE(IN  MILES) 

4 

CRITERIA  FOR  COMPUTATION  OF  DO  DEFICIT 

INPUT  VALUES 

UPSTREAM  FLOW  2.000000 

EFFLUENT  FLOW  2.000000 

UPSTREAM  BODC  3.000000 

REAREATION  RATE  COEFFICIENT  1.498322 

ENTER  VALUES  FOR  O'CONNORS  PHOTOSYNTHESIS 

ENTER  AMPLITUDE 

2 

ENTER  RATE  OF  RESPIRATION 
3 

ENTER  TAU 


2 

DISSOLVED  OXYGEN  PROFILED  FOR  DIST  OF  4.000000  MILES 


TIME  IN  DAYS 

DIST  IN  MILES 

DO  DEFICIT 

DISS  OXYGEN 

0. 000000 

0. 000000 

-5.444165 

14.285311 

0.200000 

1. 309088 

-2.965776 

11.806922 

0. 400000 

2.618176 

-1.493284 

10.  334430 
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Example  2 


! 

^ % DOSAG 

OPTION  1— USE  OBSERVED  DATA  TO  CALC  COEFF 
FOR  SINUSOIDAL  ANNUAL  WATER  TEMP  CURVE 
IF  DESIRED  TYPE  1 IF  NOT  0 
0 

OPTION  2— USE  OBSERVED  DATA  TO  CALC  CARBONACEOUS 
DEOXYGENATION  COEFF.  VIA  BOD.  TYPE  1 OR  0 
1 

OPTION  3-USE  OBSERVED  DATA  TO  CALC  NITROGENOUS 
DEOXYGENATION  COEFF.  VIA  BOD.  TYPE  1 OR  0 
0 

OPTION  4-CALC  DISSOL  OXYGEN  SAG  FOR  A SINGLE 
REACH  OF  STREAM  WITH  ONE  WASTE  DISCHARGE  VIA  DO.  TYPE  1 OR  0 
0 

OPTION  5— USE  CHURCHILL'S  METHOD  FOR  REAERATION 
OTHERWISE  USE  LANGBIEN-DURUM.  TYPE  1 OR  0 
0 

OPTION  6--INCLUDE  PHOTOSYNTHESIS  USING  CAMPS  METHOD. 

TYPE  1 IF  DESIRED,  IF  NOT  TYPE  0 
0 

OPTION  7— INCLUDE  PHOTOSYN  USING  O'CONNOR'S  METHOD 
TYPE  1 OR  0 
0 

OPTION  8--INCLUDE  NITROG  DEOXYGENATION  IN  DOSAG  CALC 
TYPE  1 OR  0 
0 

OPTION  9— COMPUTE  TEMP  GIVEN  THE  EQUATION  COEFF,  MEAN 
ANN  TEMP,  + DAY  OF  YEAR.  TYPE  1 OR  0 
0 

ENTER  TOTAL  # OF  BOD  OBSERVATIONS  FOR  CARBONACEOUS  DEMAND 
3 

ENTER  TOTAL  # OF  BOD  OBSERVATIONS  FOR  NITROGENOUS  DEMAND 
3 

ENTER  CARBONACEOUS  RATE  COEFF 

.2 

ENTER  NITROGENOUS  RATE  COEFF 
.4 

ENTER  TZERO 
3 


ENTER  DAY  OF  OBSERVATION  OF  GARB,  DEMAND-WHEN  FINISHED  TYPE  999 

4 

5 

6 

999 

ENTER  BOD  GARB  DEMAND  - WHEN  FINISHED  999 

4 

5 

6 

999 

ENTER  DAY  OF  OBSER.  OF  NITROGENOUS  DEMAND 
WHEN  FINISHED  TYPE  999 
12 

13 

14 
999 

ENTER  NITROGENOUS  DEMAND 
WHEN  FINISHED  TYPE  999 
7 
6 
5 

999 

RATE  GOEFFIGIENTS 

GARBONAGEOUS  BOD  PROGRESSION  INPUT  DATA 
TIME  IN  DAYS  BOD  SATISFIED 

4.000000  4.000000 

5 .000000  5 .000000 

6 .000000  6 .000000 
CARBONACEOUS  RATE  COEFF  = 0.012094 

GARB  ULT  BOD  168.773276 

COMPUTED  TOTAL  BOD  PROGRESSION 
TIME  IN  DAYS  TOTAL  BOD  SAT 
0.000000  0.000000 

0.500000  1.017458 

1.000000  2.028783 

1.500000  3.034011 

2.000000  4.033179 
END  OF  MODEL 


% DOSAG 


5 CONCLUSIONS  AND 
RECOMMENDATIONS 


This  report  has  described  a basic  water  quality  model 
which  can  be  used  to  quantify  impact  on  the  DO  con- 
tent of  a stream  resulting  from  inflow  from  Army 
activities.  The  model  is  limited  to  only  the  DO  aspects 
of  the  stream.  The  model  is  applicable  for  a single- 
reach, single-point-source  (as  opposed  to  area  sources) 
pollutant  inflow  situation;  in  addition,  the  flow  condi- 
tions in  a stream  must  correspond  to  a steady-tlow 
state.  This  water  quality  model  is  applicable  to  the 
type  of  stream  conditions  prevalent  in  military  installa- 
tion watersheds. 

It  is  recommended  that  Army  environmental  planners 
and  decision-makers  use  this  analytical  model  for  quan- 
tifying impact  of  new  Army  projects  and  activities  on 
dissolved  oxygen  aspects  of  water  quality. 
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